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Abstract. 

The Aharonov-Bohm effect is the prime example of a zero-field-strength 
configuration where a non-trivial vector potential acquires physical significance, a 
typical quantum mechanical effect. We consider an extension of the traditional A-B 
problem, by studying a two-dimensional medium filled with many point-like vortices. 
Systems like this might be present within a Type II superconducting layer in the 
presence of a strong magnetic field perpendicular to the layer, and have been studied 
in different limits. We construct an explicit solution for the wave function of a scalar 
particle moving within one such layer when the vortices occupy the sites of a square 
lattice and have all the same strength, equal to half of the flux quantum. From 
this construction we infer some general characteristics of the spectrum, including the 
conclusion that such a flux array produces a repulsive barrier to an incident low-energy 
charged particle, so that the penetration probability decays exponentially with distance 
from the edge. 
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1. Introduction 

In classical mechanics it is said that the vector potential has no physical meaning. Due 
to gauge invariance, only the electromagnetic field tensor has physical (measurable) 
effects. In quantum mechanics, however, the vector potential appears in gauge- 
invariant quantities that describe a new class of effects. In these cases, corresponding 
to topologically non-trivial configurations, we recognize the importance of the vector 
potential, even when the electromagnetic field vanishes everywhere in the regions 
accessible to a charged particle. 

The standard example of this class is the Aharonov-Bohm (A-B) effect [1], in which 
a magnetic field is confined to a region of space, and electrically charged particles are 
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only free to move outside this region. Although a particle cannot experience the field 
strength directly, the covariant momentum 



is affected by this configuration, that is, the vector potential carries a 'memory' of 
the presence of the magnetic field even outside the region where the field is localized. 
In this way, the particle is infiuenced by the field, through a shift in the phase of the 
wave function 



where $ is the total magnetic fiux inside the circuit (i.e. a closed path of the particle). 
This explains why the effect is called 'topological': the behavior of the particle is sensitive 
to the overall configuration of the system, even though there is no classical magnetic 
force at any point. 

The extension of the A-B problem in the presence of many localized fiuxes cannot 
be tackled exactly in general. There exists a simple argument [3] due to Aharonov 
which shows, using the Bloch theorem, that an infinite line of equispaced point-like 
fluxes would constitute an impenetrable barrier to a particle of sufficiently low energy. 
The particle would not be able to pass through such an array because it could not satisfy 
simultaneously on both sides of the barrier the Bloch periodicity conditions on its phase, 
in the light of the A-B effect. 

We are interested in exploring a possibly more realistic set-up by studying the 
propagation of a charged particle through a medium filled with point-like fiuxes. 

Experimentally, one might find a situation similar to this inside a Type II 
superconducting layer in the presence of a large magnetic field perpendicular to the 
layer. Quasiparticles in the layer would encounter numerous vortices, each containing 
a superconductor fiux quantum, and under some conditions might not penetrate the 
vortices (see, for instance, In the fractional quantum Hall effect, the strong magnetic 
field piercing a two-dimensional system is considered to be localized in fiux tubes similar 
to those in the superconducting scenario [5]. 

Configurations of this type have been addressed by several authors in recent years, 
especially in connection with the Hall problem. Some authors [6] have considered the 
approximation in which the motion of electrons is restricted to a two-dimensional lattice, 
and each plaquette is characterized by a different magnetic flux. Others have considered 
a 2-dimensional electron gas (2DEG) in a random distribution of vortices, in different 
regimes of both spatial and strength distribution (see [3 El El HDl ttH IISl UHl [H] for 
instance). 

There are many interesting aspects. Already the mathematical structure of the 
quantum mechanical problem and the existence of solutions and zero modes has 
fascinated many authors [HI [T21 [131 E] • While the scattering of particles on a single 
vortex has an exact solution [1], it becomes very complex already with two vortices 
|15j . For an arbitrary number of fluxes there is no general solution, with the possible 
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exception of Nambu's approach [2]. The existence of zero-modes has been addressed in 
[H] for the spinful case of a Pauh Hamiltonian. 

Here we want to concentrate on a spinless problem. Moreover, we will take the 
vortices to be pointlike. This sort of configuration has been addressed by Ouvry and 
coauthors in several papers [TU]. After some analytical preliminaries, they resorted to 
numerical methods to compute the density of states of electrons in a medium filled with 
pointlike fluxes. They found that for small to moderate flux strengths, the singularities 
are smothered and the 2DEG essentially sees an average effective magnetic field, which 
develops standard Landau levels, broadened by some disorder. In the strong field limit, 
the picture changes drastically, the topological nature of the problem cannot be avoided 
anymore, and one should observe a depletion of states for zero energy. 

The former observation is consistent with the analysis of [H]. In this work, Kiers and 
Weiss investigated the validity of a mean field treatment of the problem. The question 
is how accurately one can replace the fluxes with an extended magnetic field, while in 
fact the particles are never subject to a Lorentz force. They find that the condition of 
validity for this approximation is that the flux strength has to be small compared to the 
deflection of an incident particle. In our work, we consider the opposite limit of large 
flux. 

Our aim is to consider a 2-dimensional layer, punctured by magnetic fluxes, and to 
study the wave-function of a single scalar particle entering this medium. For simplicity, 
we take the vortices as point-like, so that the space available for the particles is a 
punctured plane. These are the same conditions Nambu applied in his work [2]. One 
of our aims is to considered a somewhat less general setting that allows a more explicit 
solution, and then compare our results to his. 

We are going to show that a lattice of impenetrable magnetic fluxes (vortices), 
such as the one described above, constitutes a barrier to a low-energy charged particle 
trying to pass through the medium. That is, the distribution of the vortices creates a 
conflguration whose topological constraints on the wave function are comparable to an 
effective repulsive potential. Qualitatively, there are a number of ways to see this: 

• The presence of the fluxes generates a non-zero vector potential inside the medium, 
raising the minimum energy (that is the square of the covariant momentum, eq. |2]) 
required for an electrically charged particle to exist in the medium, 

• Particles are repelled by the vortices, as their wave functions must vanish on the 
vortex sites. Therefore, the bigger the typical amplitude of the wave function in 
the flux- containing region, the bigger the energy due to the sharp spatial variation. 
This means that for low-energy states the wave function will not be able to reach 
a value appreciably different from zero in the presence of fluxes, 

• The analysis of Nambu in [2] indicates that the medium constitutes a barrier even 
from the point of view of angular momentum. In his aforementioned paper, he 
argues that the angular momentum of a particle should be greater than the magnetic 
flux present in the medium if the particle wave function is to satisfy the boundary 
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conditions. In other words, the lower angular momentum levels are missing and are 
not part of the spectrum. 

These arguments are corroborated by the aforementioned numerical simulations 
[To] showing a Lifschitz tail in the density of states at low energies for a random 
distribution of vortices. From a physical point of view, it seems quite clear that a 
charged particle approaching the medium with sufficiently low energy will be repelled, 
that is, its penetration will be exponentially damped. In the same way, if we localize 
a particle in its ground state in a region without vortices, the particle will not be able 
to escape outside that region through one containing vortices except by tunneling, and 
we should be able to construct a bound state of topological character (actually a very 
long-lived resonance), even though there is no classical force. The fact that a bound 
state can be topological in nature is somewhat new and was only implicitly suggested 
in the work of Nambu [2] . 



2. Mathematical preliminaries 

We concentrate on the case in which all the N fluxes have equal strength $ = $o/2, 
where $0 = 27r| is the quantum unit of flux. In this case it can be shown (see, for 
instance, [IE]) that the problem is invariant under time-reversal, and we can therefore 
choose the wave functions to be real. 

Indicating with {xi,yi), i = 1. . .N, the coordinates of the vortices, we can write 
the vector potential in the standard circular gauge as 

/ N N 

(A 4 V ^""^ 



^ (x - XiY + {y- Vif ' —"^{x- + iy- Vif 

N , 

Wj^tan-i'^"^^ 



X X ' 

^ '{x-Xj)^\{y-y^) 



^ \ {x - Xj) - i{y - yj^ I 



VxA =27r<l>^52^x-Xi,y-yi). (4) 

1=1 

The equation of motion for a particle in this medium is given by the Schrodinger 
equation (in units h = e = 1) 

-^(V-iA)'^ + E^ = 0, (5) 
2m 

and, in these units, integer values of $ correspond to an unobservable, quantized flux 
(in our case $ = 1/2, i.e., half a quantum of flux). 

Following Nambu's idea [2], we implement a singular gauge transformation G to 
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remove the vector potential: 

In this way, we reduce our problem to a free-field case 

- = Ei,, (7) 

with non-trivial (topological) boundary conditions on the wave functions in the region 
surrounding each vortex. 

In constructing our solutions, we must require that the wave functions vanish on 
the vortex sites 

ij{x = Xi,y = yi)=0 i = l...N , (8) 

and that they acquire the A-B phase e^™* = — 1 each time a particle completes a turn 
around a vortex. More precisely stated, in this singular gauge the effect of the vector 
potential is represented by a phase-matching condition on the wave function 

^/j (6) = -tfj {6 + 27^) (9) 

where 6 is the azimuthal angle about the vortex. 

We know from standard complex analysis that this condition implies the existence 
in the 2-dimensional plane of a cut connecting two distinguished Riemann sheets. For a 
real wave function this last condition implies that there exists at least one line exiting 
each vortex site on which the function has to vanish in order to change its sign. 



3. Construction of the solutions 



The construction of the solution for a general distribution of fluxes is not easily 
attainable (as argued in PJ). We do not need to confront these complications in order to 
show our point, and so we shall simplify the problem by taking the vortices as located 
on the vertices of a square lattice of lattice spacing L (see Fig. [T|, a case for which we 
shall be able to give an explicit solution to the problem. 

Inspired by a recent construction [T7j, we are going to estimate the minimal energy 
required for a charged particle to exist in the medium, and also to calculate the decay 
factor of particles with zero energy in the lattice. 

Before we construct the solution, it may be helpful to spend a few more words on 
our boundary conditions. Since we can take the wave function to be real, we translated 
its phase shift around each vortex with the condition that the solution has to vanish 
along one line, but we have not specifled this line. This line is not the familiar cut in a 
complex plane (which is, of course, a gauge choice). In fact, we have some freedom in 
the choice of the line along which the wave function vanishes, but this is not a gauge 
freedom in that it has a measurable effect. It would be better to say that the position 
of this line is a freedom of choice for the wave function. Therefore, in order to impose it 
as a boundary condition, we have to make this choice appropriately for the problem we 
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Figure 1. The vortices are Figure 2. The vortices are 

located on the sites of a square paired and connected by segments 

lattice. on which the wave function has 

to vanish in order to satisfy the 
topological conditions. The grey 
lines indicate the real periodicity 
of the lattice and identify the 
fundamental region over which we 
shall work. 

want to study (this consideration will be important when we consider the penetration 
of a zero-energy solution inside the medium). 

Let us consider for a moment just a pair of vortices. If we choose the line on 
which the wave function has to vanish as the ray exiting one vortex and pointing in 
the direction of the other one, we can see that the boundary conditions become that 
the function has to vanish only along the segment connecting the two fluxes; this is 
certainly a very convenient choice, compared to other solutions which would require 
the wave function to vanish on two semi-infinite lines and therefore to develop higher 
gradients. 

To construct the lowest energy solutions let us consider the vortices in pairs, 
connecting nearest neighbors with line segments along which the solution has to vanish. 
For definiteness, we connect fluxes on the horizontal direction, requiring the wave 
function to change sign when it crosses these segments (see Fig. |2]). 

Along these segments the wave function possesses odd parity. If we are interested 
in the low energy modes, this means that along the continuation of these segments, the 
function will be even and so its derivative must vanish there. To conclude our analysis 
on the boundary conditions, we notice that our system is clearly periodic. To ensure 
periodicity of the wave function, we require its derivative to vanish identically along the 
sides of each square centered on a flux (see Fig. |3]) . 
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Figure 3. Boundary conditions 
and parity of the wave function: 
the black continuous lines rep- 
resent Dirichlet boundary condi- 
tions, while the grey dashed lines 
indicate Neumann conditions. 




Figure 4. The region over which 
we construct the fundamental so- 
lution. The rest of the lattice can 
be covered starting from this ba- 
sic tile. The continuous black line 
indicates where the wave function 
must vanish (Dirichlet condition) 
and the dashed ones where its 
derivative is zero (Neumann con- 
dition). We expand the solution 
on a basis in the region I and on 
another basis in the region II and 
we impose continuity of the func- 
tion and derivative across the grey 
line. 



Bearing these considerations in mind, we now have to solve a problem with mixed 
Dirichlet and Neumann boundary conditions. We can further reduce the system under 
study and concentrate on two of the quadrants around a flux site, because the rest of 
the lattice can be covered by mirroring and flipping this unit (Fig. |4]in reference to Fig. 
§. 

In summary, we now have to solve the problem of a free particle in a rectangular box 
with sides of length 2 and 1 (in units of half of a lattice spacing). We impose Neumann 
boundary conditions everywhere, except on half of one of the long sides, where we require 
the Dirichlet boundary condition. 

This is a non-standard problem; as we are not aware of any previous study on a 
system with these boundary conditions, we shall proceed in constructing the solution 
starting from a basis compatible with the conditions. In region I of Fig. |4] we 
identify a convenient basis in the set {cosh [fc„(l -|- x)] cos(n7ry)}^p, while in region 
II we expand the solution on {cosh [Kn{l — x)] sin [{n + |)vry] }^g, with the condition 
n^vr^ ~ k'^ = {n + ^)^7r^ — = 2mE. Note that, in order to satisfy this condition, the 
knS and K^s are not necessarily real. 

By matching the wave function and its derivative across the line a; = 0, we may 
seek the values of e = 2mE for which the system admits a solution. In principle, this 
would involve the calculation of the determinant of an infinite matrix. To obtain an 
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Figure 5. We truncate the infinite-dimensional matrix to a size N and we find the 
first energy eigenvalue Sq = 2mEQ corresponding to this shortened system. This is the 
plot of N versus Eq and its fit with a polynomial in inverse powers of N up to the third 
order (higher orders do not contribute appreciably). 

approximate solution, we truncated the system to a finite size, and found tlie first energy 
eigenvalue Eq = 2mEQ as a function of the size of the matrix (see Fig. [s]). Then, we 
plotted Eo versus the order N of the matrix and performed a fit with a polynomial in 
inverse powers of N, taking the zeroth-order coefficient as the solution we would have 
got by considering the whole infinite system. 

In this way, we found the first energy eigenvalue for our solution to be: 

Eo = 2mEo = (1.0341 ± 0.0002) x ^, (10) 

that is 

Eo = (2.0682 ± 0.0002)m-^L-l (11) 

Next, we are interested in estimating the decay factor of a particle entering the 
medium with zero energy. This problem depends on the direction in which the particle 
is traveling, in that it is connected with the choice of the ray/segment over which the 
solution has to vanish. Depending on the direction of motion, the wave function may 
'choose' different configurations for these segments. 

We solve the problem for a particle moving along the x direction. That is, we 
construct a solution which exhibits periodic behavior in the y direction and real decay 
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□ 



Figure 6. Decay of the zero-energy solution moving horizontally. We require 
periodicity on the vertical axis and exponential decay in the horizontal direction. 
The continuous black line indicates where the wave function must vanish (Dirichlet 
condition) and the dashed ones where its derivative is zero (Neumann condition). 



in X (Fig. [6]). Again, we expand the wave function in appropriate bases: in region I 
and III of figure^ we use {e"''^ cos(n7r?/)}^g for right-moving and {e~"™ cos(n7ry)}^Q 

for left-moving modes. In region II we expand on < e'-^^a)'^^ sin [(n + ^)7r?/] > for 
right-moving and < e"*-""'"?^'^^ sin [(n + ^)vr?/] > for left-moving modes. 

I J n=0 

We impose matching of the wave function and its derivative across the lines x = —1 
and X = 1 and we write the damping of the solution by requiring an exponential 
suppression: 

^(x = -2,y) = e'^i:{x = 2,y) , ^(^ = -2,1/) = = 2,Z/).(12) 

We look for the values of K for which the system admits solution. 

As before, the system of equations is infinite-dimensional, so we found the lowest 
value for 7^ as a function of the order N of the matrix and performed a fit with inverse 
powers of to retain the zeroth order of the polynomial as the solution (see Fig. [?]). 

In this way, we find a decay factor for a particle moving along the horizontal 
direction: 

7^= (0.88 ±0.01) X I = (1.76 ±0.02)L-^ , (13) 

and the same K holds for a particle moving in the vertical direction because we have 
the freedom to rotate the system by 90 degrees and rearrange the segments connecting 
the vortices in the new direction. 

4. Conclusions and Outlook 



Considering a lattice of point-like magnetic vortices, we showed that the spectrum for 
a particle in such a medium is discrete, and that, for a finite lattice spacing, the lowest 
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Figure 7. We truncate the infinite-dimensional matrix to a size N and we find the 
lowest value for the decay factor K corresponding to this finite system. This is the 
plot of N versus K and its fit with a polynomial in inverse powers of N up to the third 
order (higher orders do not contribute appreciably). 

energy eigenvalue is greater than zero, by explicitly constructing the wave function 
(^0 = (2.0682 ± 0.0002)m-iL-2). 

This contrasts with what was predicted by Y. Nambu in |2]. In his paper, the 
author argues that a solution of the Schrodinger equation in our gauge would have to 
be either holomorphic, or anti-holomorphic. 

His argument goes as follows: let us switch to complex coordinates to describe the 
plane. The free particle equation now reads: 

d,d,ij = E^/j (14) 

and therefore the solution for zero energy is either analytical or anti-analytical. Nambu 
argues that, by continuity, this property should persist at higher energies as well. 
However, in the preceding section we constructed a nonzero-energy solution which clearly 
is neither holomorphic, nor anti-holomorphic, nor a linear combination of the two. 

The analyticity or anti-analyticity of the solutions is an important point of Nambu's 
construction that leads him eventually to conclude that the states with lower angular 
momentum are not admissible in the spectrum. This would imply that a particle 
entering the medium with zero energy would undergo a suppression which is not merely 
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exponential, but at least Gaussian. For that reason, we argue that our approximation 
comes closer to the true behavior, because by allowing more penetration it reduces 
uncertainty-principle energy. This statement applies even for zero energy, where 
Nambu's argument appears rigorous at first sight from (14). The loophole, we believe, 
is that for strong vortices not all of them have the wave function rotating in phase in the 
same direction [This reduces the net variation of the wave function, and clearly lowers 
the energy, which of course never can be less than zero]. In other words, at some of 
the vortices the wave function is analytic, and at some it is anti-analytic. Therefore the 
wave function as a whole is neither analytic nor anti-analytic. 

We computed the decay factor for a zero-energy particle moving along one of the 
lattice directions to be = (1.76 ± 0.02)L~^, and showed that this decay is purely 
exponential. The magnitude of this suppression depends on the direction of travel. To 
compute the decay factor in other directions it would be necessary to modify ad hoc 
the boundary condition (the positioning of the ray where the wave function vanishes). 
The condition we worked with is the one that minimizes the extension of such rays 
and therefore seems to pose the minimal constraint on the solution. Any other choice 
would have a greater impact on the shape of the wave function and would change the 
effective decay length. The directional dependence is easy to understand, because the 
coupling between charge and vortex is strong, so that the lattice length scale and the 
decay length are comparable: in the limit of vanishing lattice constant the decay length 
also vanishes. A quantitative analysis for generic directions would require a different 
formalism from the one implemented here. 

Nonetheless, we believe that the order of magnitude of the effect has been 
established, in that the lowest energy eigenvalue and the decay rate K for zero energy 
agree quite well, especially if one takes into account that the wave function still has a 
periodic variation along with the exponential decay. Such a solution is characterized 
by a real and an (orthogonal) imaginary wave-vector, equal in magnitude, to guarantee 
zero energy. It seems sensible that the real wave vector should be larger in magnitude 
than for the lowest-energy solution, because orthogonality to the imaginary vector is an 
extra constraint. The real decay rate K might be viewed as arising from an effective 
potential inside the medium. For instance, if we think about it from a WKB point of 
view we have 

1 1"^ 

K = - ^/V{x)dx . (15) 

This argument implies that the topological constraints imposed by the configuration 
of vortices act as an effective repulsive potential of order unity (in units with mass 
m = 1). This potential is clearly not constant, and in principle its precise value can be 
calculated from a detailed knowledge of K. It is more meaningful, however, to consider 
the average potential over a unit cell. As we just argued that K depends on the direction 
of travel of the particle, we see that this average effective potential is direction-dependent 
as well. 

The simple expectation is that the lowest energy eigenvalue has to be equal to 
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the average potential. However, the contribution of the relatively large orthogonal real 
vector mentioned above can make the imaginary wave number bigger than implied by 
equating the potential energy to the energy of the lowest solution. In our calculation, 
we found a good, but only approximate agreement. Different directions of travel would 
feel a different potential and, conceivably, generate a better agreement. The important 
result here is that the existence of an exponential decay, together with its magnitude, has 
been established and it can be interpreted as the effect of an average effective potential. 
Such a potential could be used to trap a particle in a region, just by surrounding that 
region with a medium of localized fluxes. Conceivably this could be a new form of 
trapping. 

The above discussion may be related to a "generalized Bloch theorem" which 
appears in various forms in the literature (for instance, see [IH], [IS])- The simplest 
version is that for a periodic potential the lowest positive energy also gives the imaginary 
wave numbers of lower-energy solutions, as if they were moving in a constant potential 
equal to the lowest positive energy. In our case, we find the imaginary wave vector at 
zero energy is bigger than this consideration would suggest. It is possible that some 
other ansatz would lower the wave vector magnitude, but for reasons discussed above 
we suspect that it still would be above the value that the naive generalized-Bloch- 
theorem would yield. This might mean that in the magnetic context there is a further 
generalization of the generalized Bloch theorem. This could be an interesting topic for 
further study. 

In |2] , Nambu argues that the proper description of the system would need to treat 
the vortices as dynamical objects themselves. Our formalism does not contemplate such 
an extension, and in the example of superconductor flux the inertia of the fluxes would 
be much greater than that of an electron. Thus the static-flux approximation makes 
physical sense. In [TU], the authors consider a random distribution of fluxes, but they are 
not interested in calculating the single-particle lowest energy level. However, it seems 
plausible to us that the order of magnitude of the decay length and the qualitative 
characteristics of the problem would not be very different from the ones found with 
our model. Our reason for saying this is that one could replace the random vortex 
distribution with a random distribution of short line segments on which the wave 
function vanishes, and this array surely would be equivalent to a repulsive potential 
of characteristic magnitude, leading to exponential, not Gaussian decay. 

The qualitative behavior we find is anyway in agreement with the analysis in flU\ . 
where it is established that for a flux strength of around $c ~ 0.3 — 0.4 a transition 
happens from a density of states in qualitative agreement with a Landau level picture to 
one characterized by a Lifschitz tail and a strong depletion of states at the bottom of the 
band, like the one we observe. In agreement with [U] as well, for $ > a mean-field 
approximation fails, and the behavior of a particle in a medium of pointlike vortices is 
completely different from the one we would observe if the particle moved through an 
extended average magnetic field. In fact, if we calculate the lowest Landau level for our 
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system, as if the particle would actually be subjected to a Lorenz force, we would find: 

pLandau ^ ^ ^ JL = = 0.25 m-^Z.-^ . (16) 



This is clearly a very different value from the one we found in our work (eq. 11), almost 
an order of magnitude smaller, showing that in our regime of large flux strength the 
mean field approximation is not valid. 

In finishing, we notice that the picture changes drastically if one introduces spin 
into the problem. In fact, as shown in many works including [13] for a Pauli Hamiltonian 
system, in certain cases particles with magnetic moment parallel to the magnetic field 
could occupy zero-modes. 
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